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When historically isolated populations meet during postglacial expansion, a mixed distribution of dis¬ 
tinct DNA lineages called contact zones is created. The gradual dissolution of the spatial genetic struc¬ 
tures in contact zones should be related to differences in pollen and seed dispersal, given no restriction 
on gene flow by e.g. reproductive isolation. We aimed to clarify effects of pollen dispersal modes on 
nuclear DNA (nrDNA) genetic structures of two codistributed species with different pollen dispersal 
modes, by analyzing nuclear microsatellites of the insect-pollinated Magnolia obovata and the wind- 
pollinated Carpinus laxiflora, which show highly consistent contact zone locations in terms of chloro¬ 
plast DNA (cpDNA). The genetic structure based on the nrDNA and that based on the cpDNA were con¬ 
cordant in M. obovata , but not so concordant in C. laxiflora. Pollen dispersal ability is higher in the wind- 
pollinated C. laxiflora than in M. obovata, resulting in the higher estimated pollen/seed migration ratio 
in C. laxiflora than in M. obovata. Therefore, the extent of postglacial lineage admixture in nrDNA was 
predominant in C. laxiflora. Our results suggested that differences in pollen dispersal ability may affect 
the nrDNA genetic structure between co-distributed species with common migration histories in the 
same area. 

Key words: admixture, contact zone, deciduous broad-leaved forest, gene flow, microsatellites, phylo- 
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Climate oscillations during the Quaternary 
affected geographical distribution of various 
plant and animal species (Hewitt 2004). During 
climatic oscillations, the ranges of organisms 
might have occasionally been divided into sever¬ 
al isolated populations. The genetic structure in 
extant populations may be the result of such his¬ 
torical changes in distribution (Hewitt 2000). 
Several organisms exhibit intraspecific genetic 
structures consistent with expected patterns from 
historical distributional changes suggested by pa¬ 
leontological studies (King & Ferris 1998, Abbott 
et al. 2000, Palme & Vendramin 2002, Petit et al. 


2002, McLachlan et al. 2005, Magri et al. 2006). 
Contact zones are defined as areas where isolated 
populations from different refugia meet during 
postglacial expansion and distinct DNA lineages 
are consequently able to mix (Harrison 1993, 
Taberlet et al. 1998). Although secondary con¬ 
tacts of isolated lineages may result in other evo¬ 
lutionary consequences, such as reinforcement 
(Song et al. 2009, Lee & Mitchell-Olds 2011, 
2013), the genetic structure of once separate pop¬ 
ulations in such contact zones are expected to 
gradually dissolve through postglacial mixing as 
long as gene flow between the lineages is not 
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strongly restricted by reproductive isolation or 
other factors. 

In plants, the rate of genetic mixing may be 
related to differences in dispersability of pollen 
and seeds. In particular, it is easy to expect that 
species with different modes of pollen dispersal 
would differ in dispersability. Based on data from 
studies using allozyme markers, Hamrick et al. 
(1990) reported that wind-pollinated plant spe¬ 
cies tend to show lower G s T values than animal- 
pollinated species. In other words, they suggested 
that immigration rates and gene flow in wind-pol¬ 
linated species are usually greater than in ani¬ 
mal-pollinated species. 

To examine the influence of different pollen 
dispersal modes on genetic structure, analyses 
using bi-parentally inherited nuclear genetic 
markers are indispensable because maternally in¬ 
herited chloroplast DNA (cpDNA) markers can 
migrate only through seeds, and thus are a reflec¬ 
tion of seed dispersal. In contrast, nuclear DNA 
(nrDNA) markers migrate through both pollen 
and seeds. Thus, it cannot be determined whether 
pollen or seed dispersal is more influential when 
analyses are based only on nrDNA or cpDNA 
markers. 

Recently, several phylogeographic studies of 
plants have used both nuclear and cytoplasmic 
DNA markers (Bai et al. 2010, Sakaguchi et al. 
2012, Qi et al. 2012, Ohtani et al. 2013, Bai et al. 
2014). For example, in wind-pollinated Juglans 
mandshurica Maxim, and ./. cathayensis Dode, 
the genetic structure observed on the basis of 
nrDNA SSR markers was not as informative as 
studies based on cpDNA markers (Bai et al. 2010, 
2014). In contrast, in insect-pollinated Kalopanax 
septemlobus (Thunb.) Koidz. and Shorea lepro- 
sula Miq., similar genetic structures were ob¬ 
served based on both nrDNA and cpDNA mark¬ 
ers (Sakaguchi et al. 2012, Ohtani et al. 2013). 
Those results suggest that differences in pollen 
dispersability between wind- and insect-pollinat¬ 
ed species are significantly responsible for differ¬ 
ences between cpDNA and nrDNA genetic struc¬ 
tures. 

The previous studies, however, did not com¬ 
pare the genetic structure between codistributed 
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species with different pollen dispersal modes in 
the same area. Therefore, the effects of dispersal 
mode on genetic structures remains unclear. We 
aimed to clarify the effects of pollen dispersal 
mode on the nrDNA genetic structure of several 
codistributed species that have different dispersal 
modes but similar distribution ranges and similar 
cpDNA genetic structures. 

Four deciduous broad-leaved tree species in 
western Honshu, Japan, Carpinus laxiflora 
(Siebold & Zucc.) Blume, C. tschonoskii Maxim., 
C. japonica Blume, and Magnolia obovata 
Thunb., with various seed and/or pollen dispersal 
modes, have been reported to show similar east- 
west genetic differentiation patterns in their cpD- 
NAs (Iwasaki et al. 2010, 2012). Furthermore, 
Tono et al. (2015) reported that the contact zones 
of the four species were highly consistent with 
each other in terms of their cpDNAs. These ob¬ 
servations indicate that the migration histories 
and extent of postglacial genetic admixture 
through seeds were similar among the four spe¬ 
cies. These species therefore provide us with a 
suitable system for examining the influence of 
different pollen dispersal modes on genetic struc¬ 
ture. 

In this study, we focused on two species, the 
insect-pollinated Magnolia obovata and the 
wind-pollinated Carpinus laxiflora , which are a 
subset of the four species examined in the previ¬ 
ous study based only on cpDNA markers (Tono et 
al. 2015). To examine the influence of different 
pollen dispersal modes on genetic structure, we 
compared the geographic patterns of genetic dif¬ 
ferentiation based on nuclear microsatellite and 
cpDNA markers between the insect-pollinated 
M. obovata and the wind-pollinated C. laxiflora. 
We addressed the following two questions: (1) 
Are the geographic patterns of genetic differen¬ 
tiation based on nrDNA markers similar to those 
based on cpDNA markers? (2) Are the geograph¬ 
ic patterns of genetic differentiation based on 
nrDNA markers different between the two spe¬ 
cies with different pollen dispersal modes? 
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Materials and Methods 

Sampling of the plants 

We used silica gel to dry leaf samples from 
250 and 264 individuals of Magnolia obovata and 
Carpinus laxiflora from 45 and 43 populations, 
respectively, in the Kinki-Chugoku region (Ap¬ 
pendix 1 and Fig. 1). Of these, 224 and 232 sam¬ 
ples from 42 and 38 populations, respectively, 
had already been examined in a previous study 
(Tono et al. 2015). Thus, 26 and 32 individuals 
from 15 and 10 populations, respectively, were 
newly sampled in this study. Voucher specimens 
have been deposited in the Makino Herbarium of 
Tokyo Metropolitan University (MAK). 

cpDNA sequencing and haplotype classification 

Twenty-six samples of Magnolia obovata and 
32 samples of Carpinus laxiflora were newly ana¬ 
lyzed in this study. In addition, the cpDNA haplo¬ 
type data from the two species were cited by 
Tono et al. (2015). The methods for DNA extrac¬ 
tion, polymerase chain reaction (PCR) amplifica¬ 
tion of cpDNA noncoding regions, sequencing, 
and classification of the eastern and western hap- 
lotypes followed those of Tono et al. (2015). The 
primers used for amplifying and sequencing the 
cpDNA non-coding regions are shown in Appen¬ 
dix 2. In total, cpDNA haplotypes were deter¬ 
mined for 250 and 264 individuals of M. obovata 
and C. laxiflora from 45 and 43 populations, re¬ 
spectively. 

Development of microsatellite markers for Carpi¬ 
nus laxiflora 

Microsatellite markers were developed 
through pyrosequencing with enriched DNA li¬ 
braries. Total DNA was extracted from the leaf 
samples of Carpinus laxiflora (voucher: AT3057) 
using the DNeasy Plant Mini Kit (QIAGEN, 
Hilden, Germany). The DNA was shotgun se¬ 
quenced using the Roche 454 GS Junior using a 
GS Junior Titanium Sequencing Kit (Roche Ap¬ 
plied Science, Penzberg, Germany) according to 
the manufacturer’s instructions on quarter plates 
with libraries of three other species identified by 


molecular identifier (MID) tags. From the ob¬ 
tained sequences, candidates of microsatellite 
loci were screened using MSATCOMMANDER 
(Faircloth 2008). The sequences with either >8 
dinucleotide repeats, >8 trinucleotide repeats, >6 
tetra nucleotide repeats, >6 pentanucleotide re¬ 
peats, or >6 hexanucleotide repeats were searched 
using Primer3 program (Rozen & Skaletsky 
2000) embedded in MSATCOMMANDER. For 
all loci, a 19-bp M13 tail (5’-CACGACGTTGTA- 
AAACGAC-3’) was added to the 5’ end as the 
forward primer sequence following the method of 
Schuelke (2000). 

PCR amplification tests were performed for 
eight individuals from six populations in a final 
volume of 5 pL using the standard protocol of the 
QIAGEN Multiplex PCR Kit (QIAGEN) for sin- 
gle-plex PCR (one primer pair per reaction). Each 
reaction contained 0.2 pM reverse primer, 0.1 pM 
FAM-labeled M13 primer, and 0.1 pM forward 
primer. PCR amplification included initial dena- 
turation at 95°C for 15 min; 30 cycles of reactions 
at 94°C for 0.5 min, 57°C for 1.5 min, and 72°C 
for 1 min; and a final extension at 60°C for 30 
min. PCR products were analyzed using an ABI 
PRISM 3100 Genetic Analyzer (Applied Biosys¬ 
tems, Foster City, California, USA). PCR product 
sizes were determined by comparison with Gen- 
eScan 500 or 600 LIZ Size Standard (Applied 
Biosystems) using GeneScan analysis (Applied 
Biosystems). The results were analyzed using 
GeneMapper version 4.0 (Applied Biosystems). 

The successfully amplified loci were ana¬ 
lyzed for 31 individuals of Carpinus laxiflora 
from one population (Gonami Pass, Ooi, Fukui 
Prefecture, Japan) with the same protocol used in 
the amplification test described above. For each 
locus, the number of alleles (AC), observed het¬ 
erozygosity ( H 0 ), expected heterozygosity (H E ), 
and fixation index (F IS ) were calculated using Ge- 
nAlEx version 6.5 (Peakall & Smouse 2006). The 
Hardy-Weinberg equilibrium (HWE) was tested 
for each locus using GenAlEx version 6.5 (Peak- 
all & Smouse 2006) with Bonferroni corrections. 
The linkage disequilibriums (LDs) between loci 
were examined using exact tests in FSTAT ver¬ 
sion 2.9.3.2 (Goudet 2002) with Bonferroni cor- 
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Fig. 1. Sampling sites in the Kinki-Chugoku region. Detailed information for each site is shown in Appendix 1. 
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Microsatellite genotyping 

For Magnolia obovata, seven microsatellite 
markers developed by Isagi et al. (1999) (M10D6, 
M6D3, M6D4, M15D5, M10D3, M6D10, and 
M17D5) were used to determine the microsatel¬ 
lite genotypes of 250 individuals from 45 popula¬ 
tions. The individuals were obtained from the 
same samples that were used in the cpDNA anal¬ 
yses. For Carpinus laxiflora, 11 microsatellite 
markers developed in the present study (Carp2, 
Carp3, Carp6, Carpi 1, Carpl3, Carpl6, Carpl7, 
Carp26, Carp27, Carp34, and Carp41) were used 
to determine the microsatellite genotypes of 264 
individuals from 43 populations. These individu¬ 
als were also obtained from the same samples 
used in the cpDNA analyses. Primer information 
is shown in Appendix 3 and 4. For each primer 
set, the forward primer was synthesized with a 
tag sequence (FAM = 5’-CACGACGTTGTA- 
AAACGAC-3’; NED = 5’-CTATAGGGCACGC- 
GTGGT-3’ VIC = 5’-TGTGGAATTGTGAGC- 
GG-3’) added to its 5’ end using the method of 
Schuelke (2000). PCR amplification was per¬ 
formed with a final volume of 5 pL using the 
standard protocol of the QIAGEN Multiplex PCR 
Kit (QIAGEN) with two- or three-plex PCR (two 
or three primer pairs per reaction). Three primers 
were used in each reaction: 0.066 pM reverse 
primer, 0.033 pM fluorescently labeled primer, 
and 0.033 pM forward primer. The methods for 
PCR amplification and determination of the prod¬ 
uct size were the same as those for microsatellite 
marker development, except for different anneal¬ 
ing temperatures (50-57°C). 

Microsatellite data analyses 

To check whether each locus met the require¬ 
ments for neutral evolution, population genetic 
analysis was performed. At first, neutral evolu¬ 
tion of the analyzed microsatellite markers was 
tested for both species using the f' sr outlier ap¬ 
proach (Beaumont & Nichols 1996, Beaumont 
2005) implemented in LOSITAN (Antao et al. 
2008). In this analysis, the relationship between 
the expected heterozygosity (H E ) and Wright’s 


5 

coefficient of genetic differentiation ( F SJ ) was ex¬ 
amined to identify outlier loci with excessively 
high or low F SJ values compared with the neutral 
expectation under an island migration model. 
This test was performed using an infinite allele 
model with 50,000 simulations, and a confidence 
interval of 0.95. LOSITAN analysis was per¬ 
formed for each population that contained more 
than 15 individuals. Second, for Magnolia obo¬ 
vata, we tested for departure from HWE at each 
locus for each population that contained more 
than nine individuals using an exact test in Ge- 
nAlEx version 6.5 (Peakall & Smouse 2006) with 
Bonferroni corrections. LD between loci was 
also tested using an exact test in FSTAT version 
2.9.3.2 (Goudet 2002). Bonferroni corrections 
were applied to all multiple statistical tests. 

To evaluate genetic differentiation in the two 
species of trees, Gst (Nei 1987) was calculated 
using FSTAT version 2.9.3.2 (Goudet 2002). In 
addition, G’ ST and standardized values of G ST 
were calculated manually (Hedrick 2005). To elu¬ 
cidate the characteristics of each population that 
contained more than four individuals, the follow¬ 
ing parameters were calculated on the basis of the 
multilocus genotype data using Gen AlEx version 
6.5 (Peakall & Smouse 2006) and FSTAT version 
2.9.3.2: allelic richness ( AR ), the total number of 
private alleles (PA), observed heterozygosity 
( H 0 ), and expected heterozygosity (Hi ). The val¬ 
ues for AR and PA were calculated by rarefying to 
five individuals using FSTAT version 2.9.3.2 and 
HP-RARE version 1.1 (Kalinowski 2005), re¬ 
spectively. 

STRUCTURE analysis 

The intraspecific genetic structure of the two 
species was estimated by Bayesian clustering us¬ 
ing STRUCTURE version 2.3.4 (Pritchard et al. 
2000, Falush et al. 2003, 2007). This analysis can 
reveal the population structure with a priori as¬ 
signment of individuals to populations on the ba¬ 
sis of multilocus genotype data. The assignment 
probabilities of clusters were estimated using the 
Markov Chain Monte Carlo (MCMC) method, 
assuming gene pools with the least possible LD 
and smallest possible departure from HWE. The 
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admixture model and LOCPRIOR model with 
correlated allele frequencies were used (Falush et 
al. 2003, Hubisz et al. 2009). Ten independent 
runs were performed using a burn-in period of 
100,000 and posterior probabilities were obtained 
from 100,000 iterations for each number of clus¬ 
ters, K - 1 to 10. The most likely value of K was 
assessed using the AK values (Evanno et al. 
2005). 


Pollen/seed migration ratios 

To estimate the relative rates of pollen and 
seed migration among populations, the pollen/ 
seed migration ratio (r) was calculated using the 
following equation: (pollen flow/seed flow) = 
{(1/G’sTn - 1) - 2 (1 /G’ STc - 1)}/(1/G’ STc - 1), where 
G’ STc is the G’ st of cpDNA and G’ STn is the G’ ST of 
nrDNA. This is a modification of the equation by 
Ennos (1994) with the substitution of G s t values 


(A) Geographic distribution pattern 
(a) Magnolia obovata 
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Fig. 2. (A) Geographic distribution pattern of two clusters estimated by STRUCTURE analysis on basis of nuclear microsatel¬ 
lite data and composition of chloroplast DNA haplotypes for Magnolia obovata (a) and Carpinus laxiflora (b). For nuclear 
microsatellite data, blue represents Cluster I; red represents Cluster II. For chloroplast DNA, eastern and western haplo¬ 
types are in blue and red, respectively. Number of samples analyzed per population and haplotype composition is propor¬ 
tional to circle size. Population numbers correspond to those in Appendix 1. 
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for F SJ values. The values of G ST and G’ sT (stan¬ 
dardized values of G sT ; Hedrick 2005) of cpDNA 
were calculated by PERMUT version 2.0 (http:// 
www.pierroton.inra.fr/genetics/labo/Software/) 
(Pons & Petit 1996) and manually, respectively. 

Results 

Geographic distribution patterns of cpDNA liap- 
lotypes 

In Magnolia obovata, five haplotypes (A, B, 
C, F, and G) were found. Three (A, C, and F) and 
two haplotypes (B and G) were classified as east¬ 
ern and western haplotypes, respectively, accord¬ 
ing to Tono et al. (2015). In Carpinus laxiflora, 
two haplotypes (A and B) were found. Haplotype 
B and A were classified as eastern and western 
haplotypes, respectively, according to Tono et al. 
(2015). All of the haplotypes observed in the two 
species had already been reported in Tono et al. 
(2015). The geographic distribution patterns of 
the eastern and western cpDNA haplotypes in the 
two species are shown in Fig. 2. After combining 
the newly obtained data from this study with data 
cited in Tono et al. (2015), proportions of the east¬ 
ern and western haplotypes still differed between 
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the eastern and western parts of the Kinki-Chu- 
goku region. The general tendencies were the 
same as in previous studies (Iwasaki et al. 2012, 
Tono et al. 2015). 

Microsatellite marker development for Carpinus 
laxiflora 

The de novo pyrosequencing produced 15,602 
reads with an average length of 435.62 bp. Three 
hundred and sixty three sequences with simple 
sequence repeats were found by MSATCOM- 
MANDER. Primer sets were successfully de¬ 
signed for a total of 41 microsatellite loci. Of 
these, 14 primer pairs (Carp2, Carp3, Carp6, 
Carpi 1, Carpl3, Carpl6, Carp 17, Carp20, Carp21 
Carp24, Carp26, Carp27, Carp34, and Carp41) 
successfully amplified DNA fragments through 
PCR (Appendix 4). The remaining 27 primer 
pairs failed to amplify DNA fragments with the 
expected product size. Characteristics of these 14 
microsatellite markers are shown in Appendix 5. 
Only one marker (Carp21) showed a significant 
deviation from HWE for the test population (P < 
0.05 after Bonferroni correction). Significant LD 
was not detected for any pairs of loci. It was dif¬ 
ficult to score fragment sizes for two markers 
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(B) Assignment probabilities 

(a) Magnolia obovata 
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Fig. 2. (B) Assignment probabilities into two STRUCTURE clusters (top) and cpDNA haplotypes (bottom) for each individu¬ 
al. Blue represents proportion of cluster I; red represents cluster II. In chloroplast DNA results, blue represents eastern 
haplotypes; red represents western haplotypes. Numbers at bottom represent population numbers. Populations are sorted 
from east to west. 
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Table 1. Parameters of genetic differentiation of chloroplasts, nuclear DNA polymorphisms, and pollen/seed migration ratios. 



Magnolia obovata 

Carpinus laxiflora 

Chloroplasts G s t (G S t c ) 

0.400 

0.675 

Chloroplasts G’ S t (G’ st t) 

0.675 

0.827 

Nuclear G ST (G st „) 

0.044 

0.029 

Nuclear G’st (G’ st „) 

0.281 

0.124 

Pollen /seed migration ratio r 

3.331 

31.608 


(Carp20 and Carp24) from 264 individuals in 43 
populations. Therefore, these three markers (Carp20, 
Carp21, and Carp24) were excluded from further 
genetic analyses. Consequently, 11 nuclear mic¬ 
rosatellite markers were used for genetic analy¬ 
ses. 

Characteristics of the microsatellite loci of Mag¬ 
nolia obovata and Carpinus laxiflora 

In Magnolia obovata, only one marker 
(M6D10) showed significant deviation from 
HWE in one (No. 40) of the eight populations ex¬ 
amined (P < 0.05 after Bonferroni correction). 
Significant LD was not detected for any pairs of 
loci in any populations. No outlier loci were de¬ 
tected with excessively high or low h'sT values 
compared with the neutral expectation. Based on 
these results, all seven loci were used for further 
analyses. 

In Magnolia obovata, AR ranged from 3.94 to 
5.5, PA from 0 to 0.24, H 0 from 0.667 to 0.833, 
and H e from 0.7 to 0.84 (Appendix 6). The mean 
values of AR, PA, H 0 , and H E were 4.87, 0.10, 0.71, 
and 0.71, respectively. 

In Carpinus laxiflora, AR ranged from 3.94 to 
5.5, PA from 0 to 0.31, H 0 from 0.6 to 0.795, and 
H e from 0.634 to 0.773 (Appendix 7). The mean 
values of AR, PA, H 0 , and H E were 5.59, 0.18, 0.76, 
and 0.78, respectively. No outlier loci were de¬ 
tected with excessively high or low T'st values 
compared with the neutral expectation. 

The overall genetic differentiation among 
populations at the 7 and 11 loci for Magnolia obo¬ 
vata and Carpinus laxiflora, respectively, was 
low (G st , G ' st = 0.044, 0.281 and 0.029, 0.124 
among the M. obovata and C. laxiflora popula¬ 
tions, respectively) (Table 1). Magnolia obovata 


had higher G ' sT values than C. laxiflora. 

Genetic structure estimated using nuclear micro¬ 
satellite markers 

In both species, the value for AK was the high¬ 
est when K = 2 (Appendix 8). Therefore, we con¬ 
sidered the appropriate number of clusters ( K ) to 
be two. The assignment probabilities into two 
clusters for each individual are shown in Fig. 4 
with information from the cpDNA haplotypes 
(Tono et al. 2015). Possible admixing individuals 
with nearly the same levels of probabilities of 
cluster I and II were observed in several popula¬ 
tions of the two species. 

Geographic distributions of the STRUC¬ 
TURE clusters in the two species are shown in 
Fig. 2, together with those of the cpDNA haplo¬ 
types. In Magnolia obovata, the proportion of 
cluster I (in blue) was high in populations from 
the eastern and Sea of Japan side of the Kinki- 
Chugoku region, whereas the proportion of clus¬ 
ter II (in red) was high in the populations from 
the western and Pacific Ocean side of the region. 
In C. laxiflora, the proportion of cluster I (in blue) 
was high in populations from the eastern part of 
the Kinki-Chugoku region, whereas the propor¬ 
tion of cluster II (in red) was high in populations 
from the western part of the region. 

Pollen/seed migration ratio 

The overall genetic differentiation among 
populations (G’ sT ) at cpDNA of Magnolia obo¬ 
vata and Carpinus laxiflora was 0.675 and 0.827, 
respectively. The pollen/seed migration ratio r of 
M. obovata and C. laxiflora estimated on the ba¬ 
sis of the values of G' ST was 3.33 and 31.60, re¬ 
spectively (Table 1). 
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Discussion 

Two nuclear genetic clusters were detected in 
both species by STRUCTURE analysis. In Mag¬ 
nolia obovata, the results based on nuclear micro¬ 
satellites and cpDNA haplotypes were highly 
consistent (Fig. 2). In contrast, in Carpinus laxi- 
flora, the genetic structure based on nuclear mic¬ 
rosatellites and cpDNA haplotypes was not as 
concordant as in M. obovata (Fig. 2). Moreover, 
in M. obovata, the difference in frequencies in 
clusters 1 and 2 was observed not only between 
the eastern and western parts but also between 
the northern (Sea of Japan side) and southern 
parts of the Kinki-Chugoku region. In previous 
phylogeographic studies of M. obovata (Iwasaki 
et al. 2012), three major cpDNA haplotypes with 
different main distribution ranges were found: 
haplotype A in eastern Japan, B in southwestern 
Japan, and C in the Sea of Japan side of Honshu. 
They concluded that the genetic structure could 
have been shaped by isolation of the populations 
into several different refugia in each of the areas 
during the last glacial period. In this study, indi¬ 
viduals with the nuclear cluster I (blue) and cpD¬ 
NA haplotype C and those with the nuclear clus¬ 
ter II (red) and cpDNA haplotype B were fre¬ 
quently observed in the Sea of Japan side area and 
western Japan, respectively. Therefore, the north- 
east-southwest genetic divergence of the nrDNA 
observed in the Kinki-Chugoku region in this 
study may be mainly attributed to historical iso¬ 
lation of the populations between the Sea of Japan 
side refugia and the southwestern Japan refugia 
during FGM. However, further studies with wid¬ 
er sampling are necessary to reveal the migration 
history that shaped the northeast-southwest dif¬ 
ferentiation in this region. 

Maternally inherited cpDNA can only be dis¬ 
persed through seeds, whereas bi-parentally in¬ 
herited nrDNA can be dispersed through seeds 
and pollen. Therefore, if geographical distribu¬ 
tion patterns for nrDNA clusters and cpDNA hap¬ 
lotypes are consistent, it means that dispersabili- 
ty through seeds and through pollen differs little. 
Additionally, the pollen/seed migration ratio in 
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Carpinus laxiflora ( r = 31.6) was much higher 
than in Mgnolia obovata (r = 3.33) (Table 1). At 
the same time, the geographic distribution pat¬ 
terns of cpDNA haplotypes in the two species re¬ 
ported by Tono et al. (2015) were very similar, 
suggesting that migration rates through seeds dif¬ 
fer little between the two species. The large dif¬ 
ferences in nuclear genetic structure as well as in 
the values of the pollen/seed migration ratios be¬ 
tween the two species likely reflect differences in 
the amount of gene flow through pollen. 

Dissolution of the genetic structure after the 
FGM may be related to differences in dispers- 
ability through pollen and seeds. Hamrick et al. 
(1990) indicated that wind-pollinated species 
have a higher rate of gene flow than do animal- 
pollinated species based on differences in G s t 
values between wind- and animal-pollinated spe¬ 
cies. For example, wind-pollinated Fagus crena- 
ta Blume was reported to show marked genetic 
structure in a molecular phylogeographic analy¬ 
sis based on cpDNA markers (Fujii et al. 2002). 
However, the genetic structure based on nrDNA 
SSR markers was not as clear as studies based on 
cpDNA markers (Hiraoka & Tomaru 2009). Sim¬ 
ilarly, in wind-pollinated Betula maximowic- 
ziana Regel, Picea jezoensis Maxim., Cercidi- 
phyllum japonicum Siebold & Zucc., Juglans 
mandshurica, and J. cathayensis, the genetic 
structures observed on the basis of nrDNA SSR 
markers were also not as clear as those based on 
cpDNA markers (Tsuda & Ide 2005, 2010, Aiza- 
wa et al. 2007, 2009, Bai et al. 2010, 2014, Qi et 
al., 2012). In contrast, insect-pollinated Melam- 
podium leucanthum Torr. & A. Gray, Sinopodo- 
phyllum hexandrum (Royle) T. S. Ying, Cyanan- 
thus delavayi Franch., Kalopanax septemlobus, 
Shorea leprosula, and Callicarpa japonica 
Thunb., showed similar and marked genetic 
structure based on both nrDNA and cpDNA 
markers (Rebernig et al. 2010, Fi et al. 2011, 
2012, Sakaguchi et al. 2012, Ohtani et al. 2013, 
Hirano et al. 2014). 

The results concordantly suggest that differ¬ 
ences between the genetic structures revealed by 
cpDNA and nrDNA markers may be caused by 
the movement of pollen, which is greater than 
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movement of seeds in wind-pollinated species. 
The wind-pollinated Carpinus laxiflora plausibly 
has a higher migration rate through pollens dis¬ 
persal than in the insect-pollinated Magnolia ob- 
ovata. Accordingly, differences in nrDNA genet¬ 
ic structure between M. obovata and C. laxiflora 
is likely due to difference in the mode of dispersal 
of pollen. 

Results from this study help to clarify the in¬ 
fluence of life history traits, such as mode of pol¬ 
len dispersal on dissolution of genetic structure. 
Hamrick el al. (1990) reported that the mode of 
pollen dispersal has a large influence on the ge¬ 
netic structure of plants. By considering phyloge¬ 
netic constraints, Duminil el al. (2007) suggested 
that differences in the mode of pollen dispersal 
do not affect genetic differentiation patterns. 
These suggestions by two previous studies are 
important but contradictory. Several reasons can 
be postulated for this contradiction. Duminil et 
al. (2007) did not consider the effects of differ¬ 
ences in the migration histories in their meta¬ 
analyses, although many phylogeographic studies 
have reported different spatial genetic structures 
and different migration histories for various spe¬ 
cies. Furthermore, their studies have included 
various geographic study systems whose topo¬ 
graphic factors and geographic scales are differ¬ 
ent. Such differences may influence migration 
histories and patterns of genetic differentiation. 
Thus, comparative studies focused on co-distrib- 
uted species with similar migration histories in 
the same geographic research system, such as 
contact zones as in the present study, will provide 
valuable suggestions about relationships between 
life historical traits and population genetic struc¬ 
ture. This study demonstrates that contact zones 
can be good systems for examining the effects of 
life history traits, such as pollen dispersability, on 
genetic structure. 
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Appendix 1. Location of the sampling sites and chloroplast DNA haplotype information for each population. 


No 

Locality 

Altitude 

(m) 

Collector 

Magnolia obovata 
Coordinates Number Haplotypes 
samples East West 

Carpinus laxiflora 

Number Haplotypes 
of 

samples East West 


1 

Japan, Shiga Pref., Yogo, Mt. Yokoyamadake 

1058 

T. Iwasaki 

35°36'N, 136°15'E 

2 

1 

1 




Tono et al. 2015, Original 

2 

Japan, Shiga Pref., Yogo, Mt. Gyoiti 

296 

A. Tono 

35°34'N, 136°11'E 

6 

6 


12 

11 

1 

Tono et al. 2015, Original 

3 

Japan, Shiga Pref., Yogo, Japan National Route 365 

637 

A. Tono 

35°42'N, 136°09'E 

6 

5 

1 

15 

15 


Tono et al. 2015, Original 

4 

Japan, Fukui Pref., Tsuruga, Kuroga forest road 

243 

A. Tono 

35°34'N, 136°03'E 




15 

15 


Tono et al. 2015, Original 

5 

Japan, Shiga Pref., Takashima, Mt. Akasaka 

539 

A. Tono 

35°27'N, 135°59'E 




8 

7 

1 

Tono et al. 2015 

6 

Japan, Fukui Pref., Kohama, Onyu Pass North 

434 

A. Tono 

35°24'N, 135°48'E 




20 

16 

4 

Tono et al. 2015, Original 

7 

Japan, Fukui Pref., Ooi, Gonami Pass 

413 

A. Tono 

35 0 2111, 135°42'E 




19 

9 

10 

Tono et al. 2015, Original 

8 

Japan, Kyoto Pref., Kyoto, Kita-ku, Kyotanba, Mt.Chourougatake 

555 

A. Tono 

35°19'N, DS^'E 




6 

2 

4 

Tono et al. 2015 

9 

Japan, Kyoto Pref., Ayabe, Mt Misen 

503 

A. Tono 

35°24'N,135°22'E 




6 

6 


Tono et al. 2015 

10 

Japan, Kyoto Pref., Yasaka, Mt. Taiko 

659 

T. Iwasaki 

35°41'N, 135°12'E 

3 

3 





Tono et al. 2015 

11 

Japan, Kyoto Pref., Fukuchiyama, Mt Shiro 

277 

A. Tono 

35°26'N, 135°08'E 




6 

6 


Tono et al. 2015 

12 

Japan, Kyoto Pref., Fukuchiyama, Mt. Ooe 

817 

T. Iwasaki 

35°27'N, BS^T 

4 

3 

1 




Tono et al. 2015, Original 

13 

Japan, Hyogo Pref., Toyooka, Toubi Pass 

379 

A. Tono 

35°25'N, 135°02'E 




8 

5 

3 

Tono et al. 2015 

14 

Japan, Hyogo Pref., Tanba, Enoki Pass 

317 

A. Tono 

SS'Tb'N, 135°01'E 

1 


1 

3 

0 

3 

Tono et al. 2015 

15 

Japan, Hyogo Pref., Toyooka, Mt.Nishitokonoosan 

609 

A. Tono 

35°25'N, 134°55'E 

16 

11 

5 

6 

5 

1 

Tono et al. 2015, Original 

16 

Japan, Hyogo Pref., Asago, Mt. Awaga 

530 

A. Tono 

35°13'N, 134°53'E 




1 

1 


Tono et al. 2015 

17 

Japan, Hyogo Pref., Kamikawa, Japan National Route 367 

445 

A. Tono 

35°07TS1,134°50'E 




1 

1 


Tono et al. 2015 

18 

Japan, Hyogo Pref., Kamikawa, Mt. Kasagata 

366 

A. Tono 

35°04'N, 134°49'E 

4 

1 

3 




Tono et al. 2015 

19 

Japan, Hyogo Pref., Asago, Touwa Pass 

542 

A. Tono 

35°19'N, 134°47'E 

2 

1 

1 

3 

1 

2 

Tono et al. 2015 

20 

Japan, Hyogo Pref., Asago, Mt. Tarugamine 

451 

A. Tono 

35°10'N, 134°46'E 

1 


1 




Tono et al. 2015 

21 

Japan, Hyogo Pref., Yabu, Mt. Myouken 

725 

A. Tono 

35°25'N, 134°39'E 

3 

1 

2 




Tono et al. 2015 

22 

Japan, Hyogo Pref., Himeji, Mt. Seppiko 

560 

A. Tono 

35°04'N, 134°38'E 




1 

1 


Tono et al. 2015 

23 

Japan, Hyogo Pref., Kami, Mt. Sobu 

786 

A. Tono 

34°13'N, 135°59 / E 

8 

4 

4 

3 

2 

1 

Tono et al. 2015 

24 

Japan, Hyogo Pref., Shiso, Kouya Pass 

702 

A. Tono 

35°11'N, 134°35'E 

3 

1 

2 

1 


1 

Tono et al. 2015 

25 

Japan, Hyogo Pref., Ichinomiya, Mt. Hujinashi 

923 

T. Iwasaki 

35°16'N, 134°35'E 

2 

1 

1 




Tono et al. 2015 

26 

Japan, Hyogo Pref., Yabu, Mt. Hyouno 

1048 

A. Tono 

35°2l'N, 134°32'E 

18 

8 

10 

6 

3 

3 

Tono et al. 2015, Original 

27 

Japan, Hyogo Pref., Shiso, Mt. Kuroo 

539 

A. Tono 

35°06'N, 134°31'E 




1 

1 


Tono et al. 2015 

28 

Japan, Hyogo Pref., Shiso, Akatani valley 

486 

A. Tono 

35°13TS1,134°30'E 




2 


2 

Tono et al. 2015 

29 

Japan, Hyogo Pref., Kami, Mt. Senno 

988 

A. Tono 

35°27'N, 134°27'E 

7 

6 

1 

3 

3 


Tono et al. 2015 

30 

Japan, Tottori Pref., Yazu, Mt. Okino 

793 

A. Tono 

35°16'N, 134°20'E 

12 

7 

5 

6 

6 


Tono et al. 2015 

31 

Japan, Okayama Pref., Bizen, Mt. Hattoji 

257 

A. Tono 

34°55'N, 134°16'E 




1 


1 

Tono et al. 2015 

32 

Japan, Tottori Pref., Yazu, Mt. Nagi 

800 

A. Tono 

35°irN, 134°11' 




8 

6 

2 

Tono et al. 2015 

33 

Japan, Tottori Pref., Tottori, Japan National Route 118 

471 

A. Tono 

35°18'N, 134°07'E 

1 

1 





Tono et al. 2015 

34 

Japan, Okayama Pref., Tsuyama, Okutusgawa valley 

427 

A. Tono 

35°13'N, 134°01'E 




1 


1 

Tono et al. 2015 

35 

Japan, Tottori Pref., Tottori, Tatsumi Pass 

561 

A. Tono 

35°19'N, 134°01'E 

3 

1 

2 

5 

5 


Tono et al. 2015 

36 

Japan, Okayama Pref., Kagamino, Mt. Izumi 

551 

A. Tono 

35°13'N, 133°59'E 

7 


7 

7 


7 

Tono et al. 2015 

37 

Japan, Okayama Pref., Kume, Mt. Tenshi 

649 

A. Tono 

34°59'N, 133°53'E 

2 


2 

6 


6 

Tono et al. 2015 

38 

Japan, Okayama Pref., Tsuyama, Mt. Genjyuji 

468 

A. Tono 

35°03'N, 133°52'E 

8 


8 




Tono et al. 2015, Original 

39 

Japan, Okayama Pref., Kagamino, Mt. Oozora 

739 

A. Tono 

35°12'N, 

4 

2 

2 

6 

1 

5 

Tono et al. 2015 

40 

Japan, Okayama Pref., Maniwa, Mt. Tsuguro 

859 

A. Tono 

34°15'N,133°47°E 

12 

7 

5 




Tono et al. 2015, Original 

41 

Japan, Okayama Pref., Maniwa, Oohira Pass 

458 

A. Tono 

35°10'N, 

5 

1 

4 




Tono et al. 2015 

42 

Japan, Okayama Pref., Kibichuo Mt. Yamato 

608 

A. Tono 

34°49'N, 133°41'E 




9 


9 

Tono et al. 2015 

43 

Japan, Tottori Pref., Kouhu, Mt. Karasuga 

1206 

T. Iwasaki 

35°2LN, 133°34'E 

4 

3 

1 




Tono et al. 2015 

44 

Japan, Okayama Pref., Shinjyo, Mt. Kenashi 

1204 

T. Iwasaki 

35°14'N, 133°30'E 

6 

2 

4 




Tono et al. 2015, Original 

45 

Japan, Okayama Pref., Takahashi Mt. Tenjin 

681 

A. Tono 

34°52'N, 133°25'E 

1 


1 

5 


5 

Tono et al. 2015 

46 

Japan, Okayama Pref., Niimi, Mt. Hanami 

1005 

A. Tono 

35°09'N, 133°24'E 




1 


1 

Tono et al. 2015 

47 

Japan, Tottori Pref., Houki, Mt. Kamakura 

594 

A. Tono 

35°15'N, \33°23^E 

5 

2 

3 

14 


14 

Tono et al. 2015 

48 

Japan, Hiroshima Pref., Jinsekikogen, Japan National Route 9 

584 

A. Tono 

34°45'N, 133°22'E 

1 


1 

8 


8 

Tono et al. 2015, Original 

49 

Japan, Okayama Pref., Takahashi, Japan National Route 313 

553 

A. Tono 

34°52'N, 

2 


2 




Tono et al. 2015 

50 

Japan, Hiroshima Pref., Shobara Mt. Dougo 

1204 

A. Tono 

35°04'N, 133°14'E 

7 


7 

5 


5 

Tono et al. 2015, Original 

51 

Japan, Hiroshima Pref., Sanwa, Mt. Hoshinoko 

719 

T. Iwasaki 

34°44'N, 133°13'E 

1 


1 




Tono et al. 2015 

52 

Japan, Hiroshima Pref., Fuchuu, Mt. Take 

378 

A. Tono 

34°38'N, 133°10'E 

9 


9 




Tono et al. 2015, Original 

53 

Japan, Hiroshima Pref., Shobara Mt. Takashiburo 

383 

A. Tono 

34°48'N, 133°08'E 




2 


2 

Tono et al. 2015, Original 

54 

Japan, Hiroshima Pref., Miyoshi, Mt. Kagenobu 

453 

A. Tono 

34°45'N, 133°05'E 




1 


1 

Tono et al. 2015, Original 

55 

Japan, Shimane Pref., Unnan Mt. Tengu 

465 

A. Tono 

35°19'N,133°04'E 

1 


1 

8 


8 

Tono et al. 2015 

56 

Japan, Shimane Pref., Okuizumo, Mt. Azuma 

727 

A. Tono 

35°05'N,133°02'E 

1 


1 




Tono et al. 2015 

57 

Japan, Hiroshima Pref., Shobara Mt. Himuro 

714 

A. Tono 

35°00'N, 132°54'E 

7 


7 

6 


6 

Tono et al. 2015 

58 

Japan, Hiroshima Pref., Takano, Mt. Ooyorogi 

926 

T. Iwasaki 

35°6'N, 132°52'E 

10 


10 




Tono et al. 2015 

59 

Japan, Hiroshima Pref., Akitakata, Mt. Oojichiyama 

539 

A. Tono 

34°39'N, 132°46'E 

12 


12 

8 


8 

Tono et al. 2015, Original 

60 

Japan, Hiroshima Pref., Higashihiroshima, Mt. Noro 

714 

A. Tono 

34°31'N,132°43 / E 

11 


11 

7 


7 

Tono et al. 2015, Original 

61 

Japan, Shimane Pref., Oota, Mt. Sanbe 

560 

A. Tono 

35°09'N, 132°38'E 




4 


4 

Tono et al. 2015, Original 

62 

Japan, Shimane Pref., Onan, Mt. Kan 

410 

A. Tono 

34°56'N, 132°31'E 

8 


8 




Tono et al. 2015, Original 

63 

Japan, Shimane Pref., Misato, Japan National Route 186 

258 

A. Tono 

35°05'N, 132°30'E 

2 

1 

1 




Tono et al. 2015 

64 

Japan, Hiroshima Pref., Hiroshima, Mt. Kaiken 

680 

A. Tono 

34°37'N, 132°29'E 

4 


4 




Tono et al. 2015 

65 

Japan, Hiroshima Pref., Geihoku, Mt. Garyuu 

1081 

T. Iwasaki 

34°4LN, 132°11'E 

3 


3 




Tono et al. 2015, Original 

66 

Japan, Yamaguchi Pref., Kano, Mt. Nagano 

1008 

T. Iwasaki 

34°16'N, Bl^'E 

15 


15 




Tono et al. 2015, Original 
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Appendix 2. Information of the chloroplast DNA noncoding regions examined in this study. 


Plant 

species 

Region 

Forward primer 

Reverse primer 

Sequesnc- Sequenced 
ing primers length (bp)l 

Reference 

Magnolia 

obovata 

trnH[G\JG)-psbA intergenic region 

ACTGCCTTGATCCACTTGGC 

CGAAGCTCCATCTACAAATGG 

Reverse 

227 

Hamilton 

1999 


trnL(\JA.A)-trnF(GAA) intergenic region 

GGTTCAAGTCCCTCTATCCC 

ATTTGAACTGGTGACACGAG 

Reverse 

358 

Taberlet et 
al. 1991 


?r«r(UGU)-f/7?L(UAA)5’exon intergenic region 

AATCGAATTCTATTCCATTAAA 

CAATCGGAAATGATTCTATC 

Forward 

356 

Iwasaki et 
al. 2012 

Carpinus 

laxiflora 

trnH{G\JG)-psbA intergenic region 

ACTGCCTTGATCCACTTGGC 

CGAAGCTCCATCTACAAATGG 

Reverse 

402 

Hamilton 

1999 


/r/iT(UGU)-t/7iL(UAA)5’exon intergenic region 

CGATTAATACTCTAAAAAGAATCTAA 

GATCTAGTGTTGGTCCAATCAGA Forward 

288 

Iwasaki et 
al. 2012 


1 Parts of PCR amplified fragments were sequenced. 


Appendix 3. Characteristics of seven microsatellite loci for Magnolia obovata. 


Locus 

Repeat motif 

Primer sequences (5'-3') 

T. (°C) 

M6D3 

(CT)„ 

Fw-ACATGGATAGTCGTTGGATA 

Rv-ACCCCACTGAAGACAAACAT 

50 

M6D4 

(CA) 2 (GA), s 

Fw-CACCGTACCCTATCAGAACC 

Rv-ATTTTCAGCATCATCAGTTG 

50 

M6D10 

(GA), 4 

Fw-AAATTGTCGTCCAACCAGTT 

Rv-AAAGCAGCAAACAGGAAGAG 

53 

M10D3 

(GA) 3s 

Fw-GTCTAGTGAGCCGCAAATGG 

Rv-GTGAACAGCTTTCTTGTGAA 

53 

M10D6 

(CT)„ 

Fw-CGACGACGAAACTACTAACA 

Rv-TTAACTTGAGGTGGAATGAC 

50 

M15D5 

(GA) 16 

Fw-GATCGTTGCTGGCTCGC 

Rv-GCCGCCTGGATTATGAA 

52 

M17D5 

(GA)„ 

Fw-TGCTGCTCGAAGTTCTGAAT 

Rv-CGTGCAGTAAATCAGGATGT 

54 


Appendix 4. Characteristics of 14 microsatellite loci for Carpinus laxiflora. 


Locus 

GeneBank 
accession no. 

Repeat motif 

Primer sequences (5'-3') 

r, (°c) 

Carp2 

LC012498 

(AT) U 

Fw-CGTCTCCATTAAAGCCAGCC 

Rv-ACTCTTGGCTCCATGTCATG 

57 

Carp3 

LC012499 

(AG), 

Fw-AGTTGCATGTGGACTGAACC 

Rv-CCTTCCATAACCGCCAAAGC 

57 

Carp6 

LC012500 

(AT),, 

Fw-CATGACACATAGAGAGGAAGC 

Rv-CTCTCCCTTCAACCCTCCAC 

57 

Carpll 

LC012501 

(AT),, 

Fw-CGCCCATTCTCTCGCATAAG 

Rv-AGCATTCCCTAGTTGTGGTG 

57 

Carpl3 

LC012502 

(AG), 3 

Fw-TGTCTCTGTTCCTCACCACC 

Rv-CCTTCCAATGTCCAATGGCG 

57 

Carpl6 

LC012503 

(AT) i3 

Fw-TTTGGTGCGGGTTTACAAGG 

Rv-GCACATCCAACCCTTACGAAG 

57 

Carpl7 

LC012504 

(AT), 

Fw-CTGCCATCAATTCCTAGGGC 

Rv-TCGAAGCTTATTTCCAGGTGC 

57 

Carp20 

LC012505 

(AG), 

Fw-CTCCGGTAGTTCTGTCCTCC 

Rv-TTCTCTCCACAAACGGACCC 

57 

Carp21 

LC012506 

(AG) 10 

Fw-TTCGGACCTGGTTCTGTCTG 

Rv-TGTCCAAACGCGCTACCAC 

57 

Carp24 

LC012507 

(AAT),u 

Fw-GGAAAGACCCACCTTGACC 

Rv-CTGTCTAAGCCTCCTGAGCG 

57 

Carp26 

LC012508 

(AGAT) 6 

Fw-AATCCAATCCCTAGGTGGCC 

Rv-AGAGTTTGAGAGTGAGGTTGC 

57 

Carp27 

LC012509 

(AAT) 15 

Fw-ACCGACCAAGAGAGTGACATC 

Rv-CGAAGGTGAGGTGACGAGTG 

57 

Carp34 

LC012510 

(AC ),2 

Fw-ACGTGGCTATCTCTGAGTCG 

Rv-AGGTATGTTCGCAATGCACG 

57 

Carp41 

LC012511 

(AT)„T(AG) I2 

Fw-TGCCTCGGTTTAAGTTCAAGC 

Rv-AGAAGAACATGGCATGCTGTAG 

57 
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Appendix 5. Results of the primer screening for 14 microsatellite Appendix 7. Population genetic parameters estimated from nucle- 


markers for one population of Carpinus laxiflora. ar microsatellite data for Carpinus laxiflora. 


Locus 

Size range 

N N a 

Ho 

He 

F is 

Population No 

N 

AR 

PA 

H 0 

H e 

Carp2 

324-354 

31 13 

0.903 

0.867 

-0.041 

2 

12 

4.91 

0.09 

0.67 

0.74 

Carp3 

328-338 

31 5 

0.323 

0.377 

0.144 

3 

15 

4.85 

0.07 

0.76 

0.74 

Caip6 

189-210 

31 14 

0.839 

0.860 

0.025 

4 

15 

4.71 

0.02 

0.72 

0.72 

Carp 11 

329-349 

31 10 

0.710 

0.815 

0.130 

5 

8 

5.12 

0.25 

0.70 

0.74 

Carp 13 

229-259 

31 13 

0.774 

0.789 

0.018 

6 

20 

5.26 

0.09 

0.79 

0.77 

Carp 16 

172-192 

31 10 

0.677 

0.807 

0.161 

7 

19 

5.11 

0.04 

0.69 

0.74 

Carp 17 

229-244 

31 10 

0.806 

0.853 

0.054 

8 

6 

4.43 

0.13 

0.74 

0.66 

Carp20 

220-244 

31 11 

0.613 

0.835 

0.266 

9 

6 

4.42 

0.00 

0.67 

0.68 

Carp21 

262-266 

31 3 

0.290 

0.530 

0.452 *** 

11 

6 

4.72 

0.00 

0.73 

0.70 

Carp24 

258-288 

31 10 

0.710 

0.807 

0.121 

13 

8 

4.74 

0.04 

0.69 

0.73 

Carp26 

358-372 

31 5 

0.710 

0.611 

-0.162 

15 

6 

5.42 

0.10 

0.80 

0.74 

Carp27 

248-287 

31 13 

0.645 

0.867 

0.256 

26 

6 

5.01 

0.09 

0.68 

0.69 

Carp34 

361-379 

31 8 

0.742 

0.758 

0.021 

30 

6 

4.95 

0.03 

0.71 

0.70 

Carp41 

198-217 

31 11 

0.871 

0.843 

-0.033 

32 

8 

4.76 

0.15 

0.77 

0.69 

N a , number of alleles; H 0 , 

observed heterozygosity; H E , expected het- 

35 

5 

5.09 

0.13 

0.65 

0.67 

erozygosity; F ls , fixation index; ***P < 0.001. 



36 

7 

4.85 

0.13 

0.64 

0.72 








37 

6 

3.94 

0.01 

0.70 

0.63 








39 

6 

4.97 

0.03 

0.73 

0.71 

Appendix 6. Population genetic parameters estimated from 

nuclear 

42 

9 

4.81 

0.07 

0.69 

0.70 

microsatellite data for Magnolia obovata. 



45 

5 

4.91 

0.00 

0.60 

0.69 

Population No 

N 

AR 

PA 

Ho 

He 

47 

14 

5.01 

0.12 

0.69 

0.73 


2 

6 

5.36 

0.21 

0.74 

0.76 

48 

8 

5.07 

0.11 

0.80 

0.71 


3 

6 

5.15 

0.00 

0.79 

0.75 

50 

5 

4.73 

0.12 

0.64 

0.64 


15 

16 

5.73 

0.22 

0.77 

0.83 

55 

8 

5.50 

0.31 

0.76 

0.76 


23 

8 

6.00 

0.16 

0.71 

0.81 

57 

6 

5.08 

0.18 

0.73 

0.69 


26 

18 

5.68 

0.24 

0.83 

0.83 

59 

8 

4.66 

0.15 

0.70 

0.69 


29 

7 

4.91 

0.22 

0.76 

0.73 














60 

7 

4.47 

0.31 

0.75 

0.66 


30 

12 

5.67 

0.09 

0.67 

0.82 








36 

7 

6.16 

0.12 

0.82 

0.81 

AR, allelic richness; PA, private allelic richness; H 

E , expected het- 


38 

8 

4.97 

0.13 

0.73 

0.72 

erozygosity; H 0 , observed heterozygosity. 




40 

12 

5.29 

0.06 

0.68 

0.77 








41 

5 

6.00 

0.16 

0.71 

0.76 








44 

6 

5.53 

0.07 

0.83 

0.74 








47 

5 

6.00 

0.27 

0.83 

0.75 








50 

7 

5.65 

0.33 

0.71 

0.79 








52 

9 

5.39 

0.13 

0.83 

0.78 








57 

7 

5.51 

0.11 

0.69 

0.77 








58 

10 

5.90 

0.19 

0.79 

0.81 








59 

12 

6.42 

0.68 

0.80 

0.84 








60 

11 

5.78 

0.10 

0.73 

0.82 








62 

8 

4.47 

0.06 

0.70 

0.70 








66 

15 

5.88 

0.24 

0.75 

0.81 








AR, allelic richness; PA, private allelic richness; H E , expected hetero¬ 
zygosity; H Q , observed heterozygosity. 
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Appendix 8. Results of Bayesian clustering (STRUCTURE, Pritchard et cil. 2000) of Magnolia obovata (a) and Carpinus laxi- 
flora (b). 


(a) Magnolia obovata 

-7900 


(b) Carpinus laxiflora 
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_ -8100 
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-11500 
2 -11600 
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£ -11800 
-11900 
-12000 
-12100 


10 


10 




K K 


The upper graphs show the mean In Pr (X|K) ± SD over 10 runs for each value of K. The lower graphs show AK values for each 
K (Evanno et al. 2005). 
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